Dynamics of delayed-coupled chaotic logistic maps: influence of network topology, 

connectivity and delay times 
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We review our recent work on the synchronization of a network of delay- coupled maps, focusing 
on the interplay of the network topology and the delay times that take into account the finite 
velocity of propagation of interactions. We assume that the elements of the network are identical 
(N logistic maps in the regime where the individual maps, without coupling, evolve in a chaotic 
orbit) and that the coupling strengths are uniform throughout the network. We show that if the 
delay times are sufficiently heterogeneous, for adequate coupling strength the network synchronizes 
in a spatially homogeneous steady-state, which is unstable for the individual maps without coupling. 
This synchronization behavior is referred to as "suppression of chaos by random delays" and is in 
contrast with the synchronization when all the interaction delay times are homogeneous, because 
with homogeneous delays the network synchronizes in a state where the elements display in-phase 
time-periodic or chaotic oscillations. We analyze the influence of the network topology considering 
four different types of networks: two regular (a ring-type and a ring-type with a central node) and 
two random (free-scale Barabasi- Albert and small- world Newman- Watts) . We find that when the 
delay times are sufficiently heterogeneous the synchronization behavior is largely independent of the 
network topology but depends on the networks connectivity, i.e., on the average number of neighbors 
per node. 
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I. INTRODUCTION 

A system composed of many nonlinear interacting 
units often forms a complex system with new emergent 
properties that are not held by the individual units. Such 
systems describe a wide variety of phenomena in biology, 
physics, and chemistry. The emergent property is usu- 
ally synchronous oscillations. Examples include the syn- 
chronized activity in pacemaker heart cells, the cicardian 
rhythms, the flashing on-and-off in unison of populations 
of fireflies, synchronized oscillations in laser arrays, in 
Josephson junction arrays etc. p], [2, El, El • 

The effect of time-delayed interactions, which arise 
from a realistic consideration of finite communication 
times, is a key issue that has received considerable atten- 
tion. The first systematic investigation of time-delayed 
coupling was done by Schuster and Wagner [5[ , who stud- 
ied two coupled phase oscillators and found multistability 
of synchronized solutions. Since then, delayed interac- 
tions have been studied in the context of linear systems 
0], phase oscillators 0, limit- cycle oscillators H, cou- 
pled maps @, [TO] , lasers [ll|, Q3 , neurons [H, 13, EH , 
etc. 

It is well-known that oscillators that interact with dif- 
ferent delay times can synchronize [l6|, [l?], EH, EH, HD, EI| , 
but the interplay of the interaction delays and the net- 
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work topology is still poorly understood. Here we review 
the main results of our recent work on the synchroniza- 
tion of chaotic logistic maps [22|, [23|, 0] focusing on the 
interplay of delays and topology. We show that a network 
of delay-coupled logistic maps can synchronize, for ade- 
quate coupling strength, in spite of the fact that the inter- 
actions among the maps have an heterogeneous distribu- 
tion of delay times. In the synchronized state the chaotic 
dynamics of the individual maps is suppressed and all 
elements of the network are in a steady-state, which is 
an unstable fixed-point of the uncoupled maps. This is 
in sharp contrast with the synchronized dynamics when 
the delays are homogeneous (instantaneous coupling and 
fixed-delay coupling), because with uniform delays the 
maps evolve in either periodic or chaotic orbits. 

We investigate the influence of the network topology 
considering four different types of networks: two regular 
(a ring-type and a ring-type with a central node) and two 
random (free-scale Barabasi- Albert [26[ and small-world 
Newman- Watts) [27j . We find that steady-state synchro- 
nization depends mainly on the average number of neigh- 
bors per node but is largely independent of the network 
architecture (i.e., the way the links are distributed among 
the nodes). This is also in contrast with the homogeneous 
delay case, because when the delay times are uniform the 
synchronization of the network depends strongly on the 
connection topology [10| . 

The rest of the paper is structured as follows. In 
Sec. [Ill we introduce the basic ingredients of our model. 
Subsequently, Sec. IIIII presents the results of the simula- 
tions. Finally, Sec. IIVI summarizes the main findings and 
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presents the conclusions from this investigation. 

II. NETWORKS AND INTERACTIONS MODEL 

We consider a network of N coupled maps: 

N 

Xi(t + 1) = (1 - e)f[xi(t)] + ^~J2 - T ^ (!) 

0% 3 = 1 

where t is a discrete time index, i is a discrete spatial 
index (2 = 1... iV) , f(x) — ax(l — x) is the logistic map, 
the matrix rj = (rjij) defines the connectivity of the net- 
work: rjij = rjji = 1 if there is a link between the zth and 
jth nodes, and zero otherwise, e is the coupling strength, 
which is uniform throughout the network, and Tij is the 
delay time in the interaction between the ith and jth 
nodes (the delay times r^- and r^ need not be equal). 
The sum in Eq. ([T]) runs over the b{ nodes which are 
coupled to the ith node (pi = Y^jVij)- The normalized 
pre-factor 1/bi means that each map receives the same 
total input from its neighbors. 

It can be noticed that the homogeneous steady-state 

xiit) = Xj(t) = x ,Vi,j,t, (2) 

where xo is a fixed point of the uncoupled map, xo = 
f(xo), is a solution of Eq. (pQ) regardless of the delay 
times and of the connectivity of the network. 

We will consider delay times that are either homoge- 
neous (jij = ro for all z, j) or heterogeneous. In the latter 
case we introduce a disorder parameter, c, that quantifies 
the degree of heterogeneity and allows varying the delays 
from a delta distribution (uniform delays) to Gaussian or 
exponential distributions. Specifically we consider 

i) Tij = ro + c£, where £ is Gaussian distributed with 
zero mean and standard deviation one. The delays are 
homogeneous (jij = ro) for c = and are Gaussian dis- 
tributed around ro for c ^ [depending on ro and c the 
distribution of delays has to be truncated to avoid nega- 
tive delays, see Fig. 1(f)]. 

ii) = To + c£, where £ is exponentially distributed, 
positive, with unit mean. The delays are homogeneous 
(jij = To) for c = and are exponentially distributed, 
decaying from ro for c ^ 0. 

To investigate the influence of the topology we con- 
sider four networks, two of them are regular and the 
other two are random. The regular ones are a ring of 
nearest-neighbor elements (NN) while in the second one 
we added a central node connected to all other nodes 
(ST). The random networks consist of a scale free net- 
work (SF) constructed according to the Barabasi- Albert 
method [26[ and, concerning the last one, we use the 
small- world (SW) topology proposed by Newman and 
Watts (13. 

In the next section we present results of simulations 
that show that the homogeneous state-state, Eq. (2), 
with xo being the nontrivial fixed point of the logistic 




0.5 1 0.5 1 5 10 

£ £ T 



FIG. 1: Xi vs. e (a) and (b); a 2 vs. e, (c) and (d). In Figs. 
1(a) and 1(c) the delays are distributed exponentially (ro = 0, 
c — 1.2, see text); the distribution is shown in Fig. 1(e). 
In Figs. 1(b) and 1(d) the delays are Gaussian distributed 
(to = 3, c = 2, see text); the distribution is shown in Fig. 
1(f). The inset in Fig. 1(d) shows with detail the transition 
to synchronization: a 2 decreases abruptly at e ~ 0.4, and is 
zero for e > 0.8. N = 500, a = 4, and p = 0.3. 



map, xo — 1 — 1/a, is a stable solution for adequate cou- 
pling strength and delay times. 



III. RESULTS 

The simulations were done choosing a random initial 
configuration, Xi(0) randomly distributed in [0,1], and 
letting the network evolve initially without coupling [in 
the time interval < t < max(r^)]. We present results 
for a = 4, corresponding to fully developed chaos of the 
individual maps, but similar results have been found for 
other values of the parameter a. 

With both, either heterogeneous or homogeneous de- 
lay times, if the coupling strength is large enough the 
network synchronizes in a spatially homogeneous state: 
Xi = Xj V z, j . Figures 1 and 2 display the transition to 
synchronization as e increases. In these figures the net- 
work has a small- world topology [27[ , but similar results 
have been found for other topologies, as discussed below. 
At each value of e, 100 iterates of an element of the net- 
work are plotted after transients. To do these bifurcation 
diagrams only the coupling strength, e, was varied; the 
network connectivity, rjij : the delay times, r^-, and the 
initial conditions, x^(0), are the same for all values of e. 

Figure 1(a) displays results for exponentially dis- 
tributed delays and Fig. 1(b) for Gaussian distributed 
delays; Fig. 2(a) for instantaneous interactions, and Fig. 
2(b) for uniform delays. It can be noticed that for small 
e the four bifurcation diagrams are very similar (we refer 
to this region as the "weak coupling region"); however, 
as e increases above ~ 0.1 the bifurcation diagrams be- 
gin showing some differences, and for large e they dif- 
fer drastically: xi is constant in Figs. 1(a) and 1(b), 
Xi = xo = 1 — 1/a, while Xi varies within [0,1] in Figs. 



FIG. 2: Xi vs. e (a) and (b); a 2 vs. e (c) and (d). In (a),(c) 
Tij = V in (b),(d) r^- = 3 V ij. All other parameters as 
in Fig. 1. 

2(a) and 2(b). 

To characterize the transition to synchronization we 
use the indicator 

a 2 = 1/N < $>*(*)- < x >] 2 > u (3) 

where < . > denotes an average over the elements of the 
network and < . > t denotes an average over time. Fig- 
ures 1(c), 1(d), 2(c) and 2(d) display a 2 vs. e correspond- 
ing to the bifurcation diagrams discussed above. It can 
be noticed that for large e there is in-phase synchroniza- 
tion in the four cases [a 2 = if and only if X{ (t) = xj (t) V 
z,j]; however, we remark that an inspection of the time- 
dependent dynamics reveals that the synchronized dy- 
namics is different: for heterogeneous delays the maps 
are in a steady-state, while for homogeneous delays the 
maps evolve either periodically or chaotically. It can also 
be observed that the four plots of a 2 vs. e are similar 
in the weak coupling regime (in Fig. 2(d) the network 
synchronizes also in a window of small e; this occurs for 
odd delays as reported in [lQ[). 

Next we analyze the influence of the network topology. 
Figure [3] displays a density plot of a 2 as a function of 
the coupling strength, e, and the mean number of neigh- 
bors per node, < b >= ^bi/N. The four panels corre- 
spond to the different networks mentioned above (SW, 
SF, ST and NN); in these cases the Gaussian distribu- 
tion of delays is the same in the four panels. Despite the 
differences when the number of neighbors is small, it can 
be observed that the synchronizability of the network is 
largely independent of the topology. 

We also show in Fig. [H for the sake of comparison two 
topologies (SW, NN) with a Gaussian delays distribu- 
tion. Analyzing Fig. 2] it is possible to conjecture that it 
represents an intermediate behaviour between the Ran- 
dom Delay Gaussian distribution and the Fixed Delay 
interaction. Regarding the case of odd delays, it looks 
like to remain a trace of the "island of synchronization" 



FIG. 3: (Color online) Random delays: synchronization re- 
gions for the four different networks considered and Gaussian 
distributed delays. The density plots represent the parameter 
a 2 as a function of e and b averaged over 10 realizations of 
the initial conditions (N = 100 and a = 4). Up: SW and SF, 
down: ST and NN. 




FIG. 4: (Color online) Random delays: synchronization re- 
gions for the SW, SF, ST and NN networks and exponential 
distributed delays. Same parameters as Fig. [3] 



present in the fixed delays interactions, but which is not 
present in the Gaussian Random Delay distribution. 

Concerning the influence of the delay distribution, in 
Fig. [5] and Figs. [6l|7lwe show the synchronization regions 
for the cases of instantaneous interactions and homoge- 
neous delays respectively. It has been reported [10] that 
there is a very important difference between the case with 
an even or odd homogeneous delay (Figs. [6j[7j). With an 
even delay it can be prove that there is a kind of "is- 
land of synchronization", it is a region for e (roughly 
speaking from e ~ 0.15 to 0.19 -weak coupling regime- 
) where the maps synchronize. In addition to this, we 
see a subcritical period doubling bifurcation in this re- 
gion. For an odd delay, this synchronization region dis- 
appears. Although this kind of delay, in the case of even 
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delays, secure these "islands of synchronization" there 
are some particularities to take into account. For exam- 
ple, in the case of the NN topology, for small number of 
neighbors the synchronization is not as good as in the 
rest of the topologies. Also in the upper limit of the 
weak coupling regime (e > 0.4) there are differences be- 
tween the irregular-networks (SW, SF) and the regular 
ones (NN, ST). 



No Delay, SW No Delay, SF 




FIG. 5: (Color online) Instantaneous interactions: synchro- 
nization regions for the four different networks considered. 
Same parameters as Fig. [3] Up: SW and SF, down: ST and 
NN. 

Considering all the showed delay distributions 
(Figs. I3|4I5I6|7() it is possible to conclude that the ran- 
dom delay distributions homogenizes the synchroniza- 
tion propensity in relation to the network topologies, 
and it also sets the value of the coupling constant e for 
the on-set of synchronization in a greater e compared 
with the other delay distributions. For instance, for ran- 
dom delays a "strong-coupling regime" starts roughly in 




FIG. 6: (Color online) Homogeneous delays (to = 5): syn- 
chronization regions for the four different networks consid- 
ered. Same parameters as Fig. [3] 




FIG. 7: (Color online) Homogeneous delays (to = 4): syn- 
chronization regions for the SW (up- left), SF (up-right), ST 
(down-left) and NN (down-right) networks. Same parameters 
as Fig. El 



e ~ 0.55, while for homogeneous interactions this regime 
starts before e ~ 0.4, what is more for instantaneous 
interactions (no delay at all) it is not clear determined 
these region, i.e. because it depends on the topology and 
the average number of neighbours. 

In further researches we plan to extend this inves- 
tigation about the special weak-coupling regime and 
the properties of synchronization related to the network 
topology and connectivity. There we have analyzed the 
properties briefly discussed here, considering five differ- 
ent kind of topologies [29[. 



IV. DISCUSSION AND CONCLUSIONS 

To summarize, we studied the dynamics of a network 
of chaotic maps with delayed interactions and showed 
that the network synchronizes in a spatially homoge- 
neous steady-state if the delay times are sufficiently het- 
erogeneous. This behavior resembles the so-called "am- 
plitude (or oscillator) death" phenomenon, which refers 
to the fact that under certain conditions the amplitude 
of delay-coupled oscillators shrinks to zero Q. It has 
been shown that that distributed delays are more sta- 
bilizing than fixed delays [19]. The stabilization of the 
fixed-point is also related to the multiple delay feedback 
method proposed by Ahlborn and Parlitz [28[, for stabi- 
lizing unstable steady states. We have recently compared 
the dynamics of a map of a network of N delayed-coupled 
maps with the dynamics of a map with N self-feedback 
delayed loops [23]. If N is sufficiently large, we found 
that the dynamics of a map of the network is similar 
to the dynamics of a map with self-feedback loops with 
the same delay times. Several delayed loops stabilize the 
fixed point; however, the distribution of delays played a 
key role: if the delays are all odd a periodic orbit (and 
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